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ABSTRACT 

Two systems of VLBI (Very Long Baseline Interferometry) at J. P. L. are now 
applicable to earth physics: an intercontinental baseline system using antennas 
of the NASA Deep Space Network (DSN), now observing at one- month intervals 
to determine UTI for spacecraft navigation; and a shorter baseline system called 
ARIES (Astronomical Radio Interferometric Earth Surveying) , to be used to mea- 
sure crustal movement in California for earthquake hazards estimation. The 
DSN system is now regularly observing between Gold stone, California and 
Madrid, Spain, determining the Earth's integrated spin rate from fringe fre- 
quency measurements . This system will soon be improved by adding the capa- 
bility to measure time delay and by extending die system to other stations of the 
DSN, making possible the determination of polar motion and of all three coordi- 
nates of the various intercontinental baselines. On the basis of experience with 
the existing DSN system, a careful study has been made to estimate the time and 
frequency requirements of both the improved intercontinental system and of 
ARIES. In this paper, such requirements for the two systems are compared 
and contrasted. The eventual requirements for the intercontinental system are 
a frequency stability of Af/f = 10 * 14 over an observing period of 24 hours, and 
a clock synchronization of 25 microseconds. The requirements on ARIES are 
less stringent in frequency, for reasons to be discussed in this paper. Over the 
shorter ARIES baselines, one must have a frequency stability of Af/f 3 • 10 ‘ M 
over 3 hours, and a clock synchronization of 25 microseconds to attain accuracy 
in each baseline coordinate of 3 cm, using an optimal observing strategy and 
bandwidth synthesis technique. 

INTRODUCTION 

A new system of Very Long Baseline Interferometry (VLBI) is being developed 
at the Jet Propulsion Laboratory (JPL) for earthquake hazards estimation. This 
system, called ARIES (Astronomical Radio Interferometric Earth Surveying), 
will be used to monitor crustal motion and regional uplift in such earthquake - 
prone areas a sourthem California and perhaps northern Mexico. The time and 
frequency requirements for ARIES are determined by two sets of parameters. 
First, one must specify the accuracy with which the crustal movement must be 
measured to give useful information concerning the earthquake mechanism, and 
then one must determine how many observations must be taken with what precision 
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and over how long a period of time in order to attain such accuracy. This paper 
outlines the logic by which time and frequency requirements are calculated for 
the system, and presents the results. 

An earthquake occurs when the accumulated stress build up in the solid crust of 
the earth, believed to be caused by underlying convective currents in the mantle, 
exceeds the elastic limit of the crust and is relieved by a sudden fracture. If 
two crustal blocks were to slide smoothly past one another, then they would yield 
freely to the force which moved them, so that no stress would accumulate and no 
earthquakes need occur. But such large blocks as we speak of here are far from 
rigid; rock considered on a scale of hundreds of kilometers is elastic and com- 
pressible. It is possible for such blocks to be sliding freely at one point of their 
contact and sticking elsewhere, while stress builds up in a complicated three 
dimensional pattern characteristic of any deformable solid. Then to measure 
the a<, emulated stress at any point, and so to estimate the likelihood and mag- 
nitude of future earthquakes, it is not sufficient to measure the relative motion 
of a few points only, or even the relative motion along the entire line of the fault, 
but measurements must be taken over whole areas tens of kilometers to either 
side of the line of contact, and integrated in a model of crustal stress and strain. 
To secure such measurements, a system must be devised which enables portable 
devices to be moved from benchmark to benchmark over wide expanses of moun- 
tain and desert. Compare Figure 1, which illustrates how two blocks may be 
expected to deform according to the elastic rebound model of the earthquake 
mechanism. 

The reality of southern California geology is even more complicated than the 
simple theory outlined above would indicate. Consider Figure 2, which shows 
the major faults, or lines of fracture, in this region. The most important single 
fracture is the San Andreas Fault, which runs southeast-northwest from the 
Gulf of California to Point Arena about 160 kilometers north of San Francisco. 

The motion along the San Andreas Fault, in its immediate vicinity, ranges from 
0 to 3 centimeters per year. Broadly speaking, this fault is the boundary between 
two interacting plates, or rigid sheets of rock into which the uppermost layer of 
the earth, the lithosphere, seems to be divided, plates which are sliding past 
one another along this line of fracture. However, the plate motion has torn the 
crust, not along one fault only, but along a whole system of parallel faults which 
divide southern California into numerous small blocks which are presumably all 
moving with respect to one another, at least as seen on a geological time scale 
(see Fig. 3). Furthermore, another system of faults exists in this region which 
runs almost due east- west, examplified by the Garlock Fault north of Goldstone, 
California, which is very difficult to account for theoretically. 

A very important symptom of an impending earthquake may be provided by the 
phenomenon of dilatancy. When the rock in a tectonically active region is 
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RELATIVE PLATE MOTION 



Figure 1. Schematic Diagram of Two Sideslipping Tectonic Plates Showing Deformation 
Due to Pressure and Fusion of Rock at the Boundary 




Figure 2. Simplified Fault Map of Southern California (Courtesy of Don L. 
Anderson, Seismological Laboratory, California Institute of Technology) 
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Figure 3. Schematic Diagram of Crustal Movements in Southern California 

(Courtesy of Don Anderson) 



subjected to a sufficient shearing stress, it develops a fine pattern of crac’.s, 
which ultimately fill with groundwa’ er. This microscopic cracking of the rock 
has two effects: the velocity of pressure waves (P- waves) thru the rock de- 
creases; and the rock increases in volume, which may be expected to induce a 
small regional uplift in the months or years prior to the earthquake which rup- 
tures the rock completely and thus relieves the strain. 

These geological phenomena determine the capability which a geodetic technique 
should have in order to be useful for earthquake hazards estimation. It should 
be capable of locating points of reference in 3 dimensions, so that regional up- 
lift as well as horizontal movement will be detected, to an accuracy of 3 centi- 
meters or better, in a coordinate system which permits results to be reproduced 
or changes measured over several decades. The VLBI technique is capable of 
measuring baseline vectors in 3 dimensions and with respect to an extragalactic 
frame of reference. It seems likely that the necessary accuracy can be attained. 
We confine ourselves here to the question: what are the time and frequency re- 
quirements to attain 3 centimeter accuracy ? 

The basic principles of VLBI geodesy are illustrated in Figure 4. The basic 
observable is the delay between the times of arrival of an electromagnetic wave 
at two antennas. If s; is the unit vector from the ith celestial source to either 
antenna, and if B is the vector baseline from one ante*—. «;o the other, and e 
is the speed of light in vacuo , then the i me delay for the ith observation is 
given by the equation 

B - Sj 


which can be rewritten in the form 


CTj = XjX + yjV + ZjZ. (2) 

where (X, Y, Z) are the baseline coordinates, and (Xj, y; , Zj ) are the coordi- 
nates cf the ith source. If the baseline coordinates (X, Y, Z) are not to be func- 
tions of time, then one must express Equation 2 and the source coordinates (Xj , 
y if Zj ) in a frame of reference rotating with the earth, allowing for the effects 
of the variation in the rotation of the earth and of polar motion. JPL now op- 
p’ ites an intercontinental interferometer between Gold tone, California, and 
Madrid, Spain, which determines changes in UT1 at intervals of approximately 
one month. At present, the Goldstone- Madrid interferometer is able to measure 
only the rate of change of time delay, which is proportional to a quantity called 
fringe frequency, but the system will soon be improved by adding the capability 
to measure time delay and by extending the system to other stations of the NASA 
Deep Space Network, making possible the determination of polar motion and of 
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Figure 4. Sketch and Flow Diagram of a Working VLBI System for Earth Physics 

Using a Portable Antenna 





all three coordinates of the various intercontinental baselines. It is planned to 
use the intercontinental interferometer at times not far removed from the dates 
of ARIES observations to determine UTl and polar motion with high precision 
(— 20 cm), and to determine a highly accurate catalog of celestial sources, so 
that the ARIES system, which will operate on relatively short baselines with a 
portable antenna, can input this information as known quantities. Thus, so far 
as the geometry is concerned, it will be possible to write Equation 2 above with 
only three unknowns, namely, ARIES baseline coordinates. 

Apart from the geometry, it is possible and necessary to solve for two other 
parameters. In general, the two clocks at the two ARIES antennas will not be 
perfectly synchronized, and they will not have exactly die same rate. One may 
include the effects of clock offset and frequency offset by rewriting Equation 2 
with two additional unknowns: 

cTj = XjX + y^Y + ZjZ + c - AT + cT • (Afff). (3) 

where AT and (Af/f) are tlw unknown clock offset and difference in clock rates, 
respectively, and T is time on a reference clock. Equation 3 is solvable with 
a minimum of five observations widely separated in azimuth. Notice that, in 
principle, VLBI can be used not only for geodesy but for clock synchronization 
in widely separated locations. 

The basic mathematics, then, sets no requirement for clock synchronization, 
since one solves for clock offset. However, there is a practical requirement. 

As illustrated in Figure 4, the noise signal received from a celestial source is 
recorded at each antenna on magnetic tape, and the tapes are cross-correlated 
by matching the streams of binary digits, called bit streams, on the tapes. (The 
signal is received at whatever frequency the receivers are tuned, heterodyned 
to generate a sine-wave signal of frequency low enough to be recorded on video 
tape, and then re-heterodyned by a technique called "fringe stopping" to allow 
approximately for the difference in Doppler effect between the two antenna loca- 
tions due to Hie rotation of Hie earth. The two resulting bit streams are then 
multiplied together. ) It is necessary that the time-tags on the two tapes be rea- 
sonably accurate if excessive time is not to be wasted searching tapes in the 
cross-correlation process. Experience with the Goldstone- Madrid interferome- 
ter and the Mark II recorder of the National Radio Astronomy Observatory sug- 
gests an a priori requirement of 25 microseconds on clock synchronization. 

Once the cross-correlation is made, it is possible to deduce the clock offset with 
a precision of about 25 nanoseconds using an instantaneous recorded bandwidth 
of 2 megahertz on a single channel, or better using a technique called bandwidth 
synthesis. 
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Although one solves for frequency offset. Equation 3 presumes that the offset 
is a constant over the observing period. Since it is desired that Equation 3 be 
over-determined and well-conditioned, it requires about 3 hours = 10 4 seconds 
to observe a sufficient number of sources— say, 12 sources with about 10 min- 
utes integration time on each, plus antenna moving time. A worst-case require- 
ment for the frequency standards can be set fay supposing that one neglects to 
include the term cT * (Af/f) in Equation 3, and that die whole term is taken up 
in a single baseline parameter— say X. In that case, in order not to exceed 3cm 
baseline errors, one would have 

(Af/f) = x/cT = 3cm/(3 - I0 10 cm/ sec > • I0 4 sec = I0'* 4 

Numerical simulations in which the effect of a constant frequency offset is in- 
cluded in the equations and in which sources are observed well distributed around 
die sky suggest that frequency variations of from 3 ' 10 14 to 8 - 10' 14 can be 
tolerated without exceeding a standard deviation of 3 cm in any baseline coordi- 
nate or in baseline length. We have adopted a value of (Af/f) = 3 • 10“ 14 as a 
reasonable requirement for the ARIES system. 

The requirement for the intercontinental interferometer is approximately (Af/f) - 
10 14 , and is more stringent than for ARIES chiefly because a longer observing 
period is needed for a good solution, about 10 hours. 

An attempt to compare the simpie theory outlined above with experimental data 
illustrates both general similarities and striking differences. Two interfero- 
metry experiments were conducted on two separate days over a short (16 km) 
baseline, using a 24 kHz instantaneous recorded bandwidth, between die Mars 
(64 meter) and Echo (26 meter) antennas at the Goldstone DSN complex, as part 
of a series of tests of die ARIES technique. On the first day (18 October 1972), 
x>th stations wei ■ equipped with hydrogen masers; on die second (21 November 
1972), one of die hydrogen masers was replaced by a rubidium oscillator. Each 
experiment consisted of 7 hours of observation, over which a single solution was 
made for AT and for Af/f via Equation 3. Residuals were formed of the observed 
fringe frequencies and time delays minus those calculated from the solution, and 
are displayed in Figures 5 and 6. The fringe frequency jitter for the rubidium- 
equipped station is about 0.5 millihertz, corresponding to Af/f = 2 - 10 , since 

all observations were at S-band (~ 2.3 gigahertz). This figure equals or ex- 
ceeds the expected performance of the rubidium oscillator, la the other hand, 
the jitter for the hydrogen maser equipped station is 0. 1? millihertz, corre- 
sponding to Af/f = 5 • 10 14 , which fails to meet the expected short-term stability 
of die hydrogen maser by about a factor of 5. However, in this case, the ob- 
served fringe frequency reflects other sources of error, not only oscillator 
instability; and we suggest that the dominant error source is produced by the 
ionosphere. Also, in the case of the time delay residuals, the dominant error 
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Figure 5. Residuals in Time for Two Experiments, One Using Hydrogen Masers, 

and the Other, a Rubidium Oscillator 
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Figure 6. Residuals in Frequency (Same Experiments as Figure 5) 




source was not the frequency standard, but rather the system noise (Fig. 5). 
Although the stability of the hydrogen maser exceeds that of the rubidium by 
about a factor of 20, accuracy of the solution for baseline coordinates increases 
only by a factor of 2, from 8 centimeters to 4 centimeters in baseline length, 
and the scatter of the delay residuals hardly improves at all. Nevertheless, the 
accuracy of the H-maser equipped stations is close to that required for geophys- 
ical purposes. Future improvements to ARIES will include local calibrations 
for water vapor in the troposphere and for charged particle content in the iono- 
sphere, and use of a wider recording bandwidth. 

ARIES has two advantages which the intercontinental system does nat enjoy: 

1. Since UT1, polar motion, and source locations will be input into ARIES 
as data, the list of parameters to be solved for is shorter, and the num- 
ber of necessary observations smaller. 

2. Since ARIES will operate on fairly short baselines, a much larger area 
of the celestial sphere will be mutually visible to both antennas, and a 
source list may be chosen much better distributed over the celestial 
sphere to optimize the solution for the unknowns. 


SUMMARY 

A simple theory of the time and frequency requirements for VLB I geodesy sug- 
gests that, with current equipment, a frequency stability of 3 • 10 ' ,4 and syn- 
chronization requirement of 25 microseconds will permit the determination of 
baseline coordinates to about 3 centimeter precision with 3 hours observing. 

Once a solution is obtained, the VLBI technique itself would permit clock offsets 
to be calculated with an accuracy of 25 nanoseconds at widely separated sites, 
or potentially to about 1 nanosecond using a bandwidth synthesis technique. Ac- 
tual demonstrations showed that a baseline accuracy of 4 centimeters was at- 
tained with a prototype system over a 16 km baseline utilizing H-maser frequency 
standards. 
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QUESTION AND ANSWER PERIOD 


DR. KLEPCZYNSKI: 

Are there any questions or comments? Yes Dr. Johnston from NRL. 

DR. JOHNSTON: 

I would like to make a comment on the new H. P. cesium standard. 

In cooperation with the Naval Observatory, on a recent very long baseline 
experiment, we attempted to use two cesium standards, one at each station. 

We did have some small problems with time synchronization, but we failed to 
get fringes with the cesium standards, which could possibly indicate short term 
stability on the order of one second, maybe as low as a part in tenth to the 
tenth, which is very poor. 

If anyone has any further comments on that, I would like to hear them. 

DR. KLEPCZYNSKI: 

Yes. 


I have a separate question. 

DR. KLEPCZYNSKI: 

Before we go to this question, are there any comments on the last one? Yes. 
DR. TOM CLARK: 

We have used one of the H. P. Cesium standards at Goldstone in an experiment, 
oh, several months ago, when the hydrogen maser at Goldstone was not in the 
best of health, and although we had fringes at 7. 8 and 15 gigahertz, the phase 
stability is probably worse on short terms, 10 seconds or so, than even a 
rubidium standard. Is there a possibility that some internal loop time constants 
can be changed in the cesium standards to assist us ? 

I know that Clark Wardrip here at Goddard has run considerable tests on the 
stability of the cesium. Perhaps you would want to talk to him about that, too. 
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DR. KLEPCZYNSKI: 


Is there a comment over here ? 

DR. HELWIG: 

At the National Bureau of Standards, we have some results with a Hewlett- 
Packard Super Tube and I think it performs according to specs, and that is all I 
can say. I think a part in 10 to the 11th is at least an order of magnitude for a 
one second averaging time. 

So this would be my comment, if you use the Super Tube just as a black box, and 
you assume that specs of H. P. are correct, I think you are in good shaoe. If 
you expect more than the specs say, you may not get it. 

DR. KLEPCZYNSKI: 

A comment from Dr. Winkler. 

DR. WINKLER: 

We will hear much more about it by Mr. Percival, who will report tomorrow on 
experiments which we have done in the observatory' with 11 units. 

There is one point in regard to Dr. Helwig's comment with which I completely 
agree. It is that I think one has to remember that the characteristics of the 
sigma, tau plot of the cesium standards are vastly different from that of a 
hydrogen maser, and of that of a rubidium standard. 

These are three different kinds of animals. The cesium standard, in general, 
will follow a one over square root of tau performance, beginning from the time 
constant of the servo loop, which is in the order of a second to a minute or 
slightly shorter than a minute down to, and that is a point of some contention, 
to the flicker level. 

It is here where the standards, according to our experience, seem to be better 
than the specifications. It is how far they will go down to long integration 
periods which determines their main quality for time keeping. 

In the short time range, if you wan* to compare them with a rubidium standard 
for periods of 100 seconds or so, they will be slightly inferior, but again I 
completely agree with Dr. Helwig, there seems to be something wrong with that 
one standard referred to before which, incidentally, had completely stopped on 
its way from the observatory, and which has not performed according to what it 
should be. 
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But more about it tori. - r row. 


DR. KLEPCZYNSKI: 

Dr. Helwig again. 

DR. HELWIG: 

As an example, we have reliable data on a comparison of one rubidium standard 
against our primary standard, we get 8 parts in 10 to the 12th one second, and 
it averages down as the square of tau, reaching a flicker level of about 2 parts 
in 10 to the 14th. 

DR. KLEPCZYNSKI: 

There was another comment over here? Yes. 

When you ask a question, please identify yourself by name and place, so that 
people will get to be familiar with everybody in the audience. 

DR. REINHARDT: 

Dr. Reinhardt, from Harvard. 

If you were to use an artificial signal source, let's say on a satellite, what 
performance could you get using hydrogen masers, and do you need that kind of 
performance? Would you benefit from that kind of performance ? 

DR. i' HEGEL: 

I imagine that we could. We have no experience using artificial sources for the 
AIRES project, of course. We had to use natural sources, for obvious reasons, 
because we want our observations finally to be reduced to an inertial system, 
which the extra galactic framework provides us, or some kind of a fixed body 
reference. 

I am certauly not ruling out the idea of using artificial sources, but we haven't 
used them up to now. 

DR. REINHARDT: 

What kind of performance did you get, signal to noise ? 
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DR. FUEGEL: 


Incredibly good, but I couldn’t put the number up. 

DR. REINHARDT: 

Could you use that kind of accuracy? 

DR. FUEGEL: 

I may have misled you on one respect. We are not really limited by the signal 
to noise ratio of the faint sources, because there are enough brighter sources to 
use to really control the solution. 

So, yes, I imagine the solution would be somewhat better, but not an order of 
magnitude better, and we are not primarily limited by the sources. It only 
appeared so because of the way that slide was displayed. 

QUESTION: 

I was just going to suggest that you look into using artificial sources, because 
these sources are all time variable, since they are very small extragalactic 
sources and they have only been studied for the last — well, they weren't known 
15 years ago, and the characteristics of these sources are still under study, so 
that perhaps using artificial satellites you would have a very good signal to noise 
which may improve your solution somewhat. 

DR. FLIEGEL: 

Yes, that is a good point. We are sticking with the quasars and the Seyfert 
galaxies primarily because we hope that we will not have the same problems as 
people who work, for example, with water vapor sources. 

We hope, for example, that the proper motion will turn out to be zero, but we 
will have to keep an eye on it. 

QUESTION: 

Can I just comment about the use of artificial satellites. We have used them and 
got results from them, but I think in connection with the particular project you 
are talking about that you do have to be concerned with the satellite orbit itself, 
as well as the librational motion of the satellite. This introduces many new 
parameters into the solution, and just from the point of view of the post process- 
ing analysis, as we call it, the analytical framework is going to be very much 
more complicated. 
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The fact, for instance, that the satellite is at a finite distance from the stations, 
and you don't have the vectors from each station parallel to each other, specifi- 
cally to the accuracy that you are talking about, is also going to have to be taken 
into account. It is a rather messy analytical problem. 

DR. FLIEGEL: 

I quite thoroughly agree. We would only use artificial sources as a last resort 
for that very reason, and they would have to be distant space probes, inter- 
planetary probes, I think rather than earth satellites. But we can talk about 
this. 

DR. KLEPCZYNSKI: 

lam going to stop the discussion on this so we can move on to our next paper, but 
before I do that I want to thank Dr. Fliegelfor v very interesting presentation. 
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